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ABSTRACT 

We present Spitzer IRS mid-infrared (5-35 fim) spectra of a complete flux- 
limited sample (> 3 mJy at 8 fim) of young stellar object (YSO) candidates 
selected on the basis of their infrared colors in the Serpens Molecular Cloud. 
Spectra of 147 sources are presented and classified. Background stars (with slope 
consistent with a reddened stellar spectrum and sihcate features in absorption), 
galaxies (with redshifted PAH features) and a planetary nebula (with high ion- 
ization lines) amount to 22% of contamination in this sample, leaving 115 true 
YSOs. Sources with rising spectra and ice absorption features, classified as em- 
bedded Stage I protostars, amount to 18% of the sample. The remaining 82% 
(94) of the disk sources are analyzed in terms of spectral energy distribution 
shapes, PAHs and silicate features. The presence, strength and shape of these 
silicate features are used to infer disk properties for these systems. About 8% of 
the disks have 30/13 /xm flux ratios consistent with cold disks with inner holes or 
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gaps, and 3% of the disks show PAH emission. Comparison with models indicates 
that dust grains in the surface of these disks have sizes of at least a few /im. The 
20 yLim silicate feature is sometimes seen in absence of the 10 /im feature, which 
may be indicative of very small holes in these disks. No significant difference is 
found in the distribution of silicate feature shapes and strengths between sources 
in clusters and in the field. Moreover, the results in Serpens are compared with 
other well-studied samples: the c2d IRS sample distributed over 5 clouds and a 
large sample of disks in the Taurus star-forming region. The remarkably similar 
distributions of silicate feature characteristics in samples with different environ- 
ment and median ages - if significant - imply that the dust population in the disk 
surface results from an equilibrium between dust growth and destructive collision 
processes that are maintained over a few million years for any YSO population 
irrespective of environment. 

Subject headings: stars: pre-main sequence - infrared: stars - stars: circum- 
stellar matter - stars: planetary systems: protoplanetary disks 



Introduction 



Newly formed stars are obs erved to have infrared (IR) exce ss due to their circumstellar 



disk composed of dust and gas (IStrom Ill992l : iHillenbrand 



20081 ). Most older main-sequence 



(MS) stars, on the other hand, have photospheric emission with no excess in the IR. It is 
intuitive to conclude that the circumstellar disk evolves with time, gradually getting rid of 
the IR excess. One of the main questions in stellar astrophysics is how this happens. 

Observational studies, as well as theoretical simulations, have demonstrated the in- 
teraction between star and disk. The stellar radiation facilitates disk evolution in terms 



1995 



of photoevaporation (e.g ., iRichling fc Yorkd l2000t [Alexander et al. 



Gorti fc Hollenbachll2009l) or dust growth and settling (e . g..lWeidenschilling 



2006: 



Alexander 



1980 



2008 



Sterzik et al 



Dominik fc TielensI 119971 : iDuUemond fc DominikI 120051 : iJohansen et al.l l2008l ) . In the 
other direction, mass is accreted from the disk to the star foll owing magnetic field lines (e.g.. 



MuzeroUe et al.ll2003l : IWhite fc Basrill2003l : iNatta et al.ll2004f ). A diversity of stellar temper- 
atures, luminosities and masses among young stars has been known and studied for decades. 
Facihtated by new IR and (sub-)millimeter observations, a great variety of disk shapes, struc- 
tures and masses is now being actively studied. The next step is to try to connect stellar 
and disk characteristics in order to understand the evolution of these systems. 



The study of a single object, however, is unlikely to provide unambiguous information 
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regarding the evolutionary stage of the associated disk. Most studies to date refer to samples 
of young stars scattered across the sky, or to sources distributed across large star-forming 
clouds like Taurus. In addition to evolutionary stage, the specific environment in which the 
stars are formed may influence the evolution of disks by dynamical and radiative interaction 
with other stars or through the initial co nditions of the starting c loud, making it difficult 
to separate the evolutionary effects (e.g., iRichling fc Yorkd Il998l . 120001 ). For this reason, 
clust ers of stars are very often used as labora tories for calibrating the evolutionary sequence 



Lada fc Ladalll995l : iHaisch et al.ll200ll ). The power of this method, to gain statistical 



(e.g. 

information on disk composition in coeval samples, was found to be very su ccessful for loose 



assoc iations of older, pre-main sequence stars such as the 8 Myr old rj Cha (iBouwman et al. 



20061 ) and the 10 Myr old TW Hydrae association ( lUchida et al.ll2004l ). Identifying clusters 
of even younger disk populations is a natural step towards the completion of the empirical 
calibration of the evolution of disks surrounding young low-mass stars. This paper analyzes 
the inner disk properties of a flux-limited, complete unbias ed sample of young stars with IR 
excess in the Serpens Molecular Cloud (d = 259 ± 37 pc, IStraizys et al.lll996l ) which has a 
mean age ~ 5 Myr, (jOliveira et al.l 120091 ) with an YSO population in clusters and also in 
isolation. It has been recently argued (L. Loinard, private communication) that the distance 
to Serpens could be considerably higher than previously calculated. This would imply a 
rather younger median age for this cloud. 

The Spitzer Legacy Program "From Molecular Cores to Planet- Forming Disks" (c2d) 
has uncovered hundreds of objects with IR excess in five star-forming clouds (Cha II, Lu- 
pus, Ophiuchus, Pe rseus and Serpens), and allowed statistical studies within a given cloud 



(lEvans et a 



Fazio et al. 



20091). The c2d study of Serp ens with IRAC ( 3.6, 4.5, 5.8 and 8.0 /xm 



20041 ) and MIPS (24 and 70 /xm, iRieke et al.l 120041 ) data has revealed a rich 



population of mostly previously unknown y oung stellar objects (YSOs ) associated with IR 
excess, yielding a diversity of disk SEDs ( jHarvey et al.l l2006l l2007al Jbl). Because of the 
compact area in Serpens mapped by Spitzer (0.89 deg^), this impressive diversity of disks 
presents itself as an excellent laboratory for studies of early stellar evolution and planet 
formation. Indeed, the Serpens core (Cluster A, located i n the northeastern part of the area 



studied by c2d) has been well studied in this sense (e.g., Zhang et al 



19921 : iTesti fc Sargentlll998l : iKaas et al.ll2004J : lEiroa et al. 



2005 



1988 



Eiroa fc Casah 



Winston et al 



2007 



2009|) 



whereas only some o f the objects in Gro up C (formerly known as Cluster C) were studied 
with ISOCAM data JPiupvik et al.ll2006l ). 



Because of its wavelength coverage, sensitivity and mapping capabilities, the Spitzer 
Space Telescope has offered an opportunity to study many of these systems (star+disk) in 
unprecedented detail. Spitzer's photometry in the mid-IR, where the radiation reprocessed 
by the dust is dominant, gives information on the shape of the disks and, indirectly, its 
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evolutionary stage (assuming an evolution from flared to flat disks). Follo w-up mid-IR 



spect roscopic observations with the InfraRed Spectrograph (IRS, 5-38, /im iHouck et al. 



20041 ) on-board Spitzer probe the physical and chemical processes affecting the hot dust in the 



provide information on dust grain size distribution and structure (e.g. 


. van . 


Boekel et al. 


2003; 


Przveodda et al. 


2003 




Bouwman et al. 


2008; 


Kessler-Silacci et al. 


2006, 


2007 




Geers et al. 


2006; 


Watson et al. 


2009; 


Olofsson et al. 


2009 


). These, in turn, reflect dynamical processes 



such as radial and vertical mixing, and physical processes such as annealing. A smooth 
strong Si-0 stretching mode feature centered at 9.8 /im is indicative of small amorphous 
silicates (like those found in the ISM) while a structured weaker and broader feature reveals 
bigger grains or the presence of crystalline silicates. Polycyclic aromatic hydrocarbon (PAH) 
features are a probe of the UV radiation incident on the di sk, whereas their abundance plays a 



crucial role in models of disk heating and chemistry (e.g.. lGeers et al.ll2006l ; iDuUemond et al. 
2007l : IVisser et al.lboOTh . The shape and slope of the mid-infrared excess provides information 
on the flaring geometry of the disks ( iDuUemond fc Dominikll2004l ). while ice bands may form 
for highly inclined sources (edge-on) where the light from the central object p asses through 
the dusty material in the outer parts of the disk (jPontoppidan et al.l 120051 ) . Thus, the 
wavelength range probed by the IRS spectra enables analysis of the geometry of individual 
disks. It also probes the temperature and dust size distributions as well as crystallinity of 
dust in the disk surface at radii of 0.1 - few AU. Statistical results from a number of sources 
help the understanding of the progression of disk clearing and possibly planet formation. 

Our group has been conducting mult i- wavelength observing campaigns of Serpens. Op- 
tical and near-IR wavelength data, where the stell ar radiation domina tes, are being used 
to characterize the central sources of these systems (jOliveira et al.ll2009l ). Effective temper- 
atures, luminosities, extinctions, mass accretion rates, as well as relative ages and masses 
are being determined. In this paper, we present a complete flux-hmited set of Spitzer IRS 
spectra for this previously unexplored young stellar population in Serpens. We analyze these 
spectra in terms of common and individual characteristics and compare the results to those 
of Taurus, one of the best studied molecular clouds to date and dominated by isolated star 
formation. A subsequent paper will deal with the full SED fitting for the disk sources in 
Serpens and their detailed analyzes. Our ultimate goal is to statistically trace the evolution 
of young low-mass stars by means of the spectroscopic signatures of disk evolution, discussed 
above. 

Section [2] describes our Spitzer IRS data. In § E] we divide our sample into categories 
based on their IRS spectra: the background contaminants are presented in § 13.11 separated 
according to the nature of the objects; embedded sources are presented in § 13.21 and disk 
sources in § 13.31 (with emphasis on PAH and silicate emission). In § H] we discuss disk 
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properties in relation to environment and to another cloud, Taurus, for comparison. In § [5] 
we present our conclusions. 



2. Spitzer IRS Data 



2.1. Sample Selection 



Harvey et al.l (j2007bl ) describe the selection criteria used by the c2d team to identify 
YSO candidates based on the Spitzer IRAC and MIPS data, band-merged with the 2MASS 
catalog. They used a number of color-color and color-magnitude diagrams to separate YSOs 
(both embedded YSOs and young stars with disks) from other types of sources, such as 
background galaxies and stars. In this manner, a set of 235 YSOs was identified in Serpens, 
in both clusters and in isolation, as described in § 14.11 These criteria, however, are not 



fail-proof. iHarvey et al.l (j2007bl ) treated this problem very carefully in a statistical sense, 
but the locus region for YSOs and, for instance, background galaxies overlap somewhat in 
color-color diagrams with Spitzer photometry. Therefore, some contamination may well still 
be present in the sample which was impossible to disentangle from photometry alone. 

An additional criterion was applied to this sample in order to guarantee IRS spectra 
with sufficient quality to allow comparative studies of solid-state features, namely S/N > 
30 on the continuum. A lower limit flux cutoff of 3 mJy at 8;Um was imposed. This flux 
limit ensures coverage down to the b rown dwarf limit (L ~ 0.01 Lq) and leads to a flnal 
sample of 147 objects, the same as in lOliveira et al.l (120091 ). This is a complete flux- limited 
sample of IR excess sources in the c2d mapped area, except in the Serpens Core. It is also 
important to note that this sample, by deflnition, does not include young stars without IR 
excess (Class III). Figure [1] shows the distribution of our observed objects in the Serpens 
Molecular Cloud. Table [T] hsts the sample together with their 2-24 /xm spectral slope. 



2.2. Observations and Data Reduction 



Out of the 147 objects in this sample, 22 were observed between September 2004 and 
June 2006, as part of the c2d IRS 2nd Look campaign. The reduced high resolution spectra 
(R = A/AA = 600; Short High [SH], 9.9 - 19.6 /im and Long High [LH], 18.7 - 37.2 /im), 
complemented with Short Low spectra (SL, R = 50 - 100; 5.2 - 14.5 /im) are public to 



the community through the c2d delivery website (lEvans et al.l 120071 ). The other 125 targets 
were observed as part of a Spitzer Space Telescope's cycle 3 program (G03 #30223, PI: 
Pontoppidan) in the low resolution module (SL and Long Low [LL], 14.0 - 38.0 /im) in 35 
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Fig. 1. — Observed objects over-plotted on the c2d IRAC4 (8.0 fim) map of Serpens. The 
contou rs (5, 10, 15, 20 an d 25 mag) of visual extinction are derived from the c2d extinction 
maps (lEvans et al.l 120071 ). Blue squares are background stars, cyan ellipses are redshifted 
galaxies, red circles are embedded sources, and stars represent disks (yellow disks with PAH 
emission and green other disk sources). The pink triangle is the PN candidate, see § IA.3l for 
details. 
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AORs, between October 2006 and April 2007. 

All objects were observed in IRS staring mode and extracted from the SSC pipeline 
(version S12.0.2 for the c2d IRS 2nd look program, and version S1 5.3.0 for the G03 pr ogram) 
basic calibration data (BCD) using the c2d reduction pipeline ( iLahuis et al.l 120061 ) . This 
process includes bad-pixel correction, optimal PSF aperture extraction, defringing and order 
matching. Technical information about each object can be found in Table [H 



3. Results 

The large number of objects studied here allows statistical studies of each stage of disk 
evolution in a single star-forming region. By simply looking at the spectra, a diversity of 
spectral shapes can be noted immediately. This diversity of objects can be divided into 
several categories, each of which is thought to be related to a different evolutionary stage of 
the star+disk system, described below. 



3.1. Background Sources 

Due to the low galactic latitude of Serpens (/=030.4733, 6=+05.1018) the contamination 
by galactic background sources is expected to be higher for this region than for other studied 
star- forming regions (e.g., Taurus, Lupus, Chamaeleon). With the help of the IRS spectra, 
background contaminants can readily be identified, even when seen through a large column 
density of molecular cloud material. These are presented and discussed in Appendix \M 

To further study the YSOs, the 32 contaminants are removed from the sample, which 
then consists of 115 objects for further analysis. 



3.2. Embedded Sources 



A newly formed protostar, still embedded in an envelope of gas and dust, will show 
an excess in the mid-IR that rises with wavelength. In this cold environment, a significant 
fraction of the molecules are found in ices. Many of these i ces are observed as absorption 
features in the wavelength range of the IRS instrument (e.g.. 



Boogert et all 120041 ) 



Of our sample, 21 sources have spectra consistent with embedded sources. These spectra 
are shown in Figure [2l The most prominent ice feature, CO2 at 15.2 /im, can readily be seen 
in all spectra. Weaker features, such as H2O at 6.0 and 6.85 /xm and CH3OH at 9.7 fim, 
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need high S/N. Those features are detected in objects #26, 27, 28, 44, 45, 47, 49, 67, 72, 73, 
135 and 140. As can be seen in Figured], the embedded objects seem to be concentrated in 
high extinction regions {Ay ~ 10 — 15 mag). 



The embedded objects observed as part of the c2d IRS 2nd Look Program have been 
st udied in detail in a series of papers on ices aro und YSOs (H2O inlBoogert et aLll2008l : CO2 

CH4 at 7.7 yum in lOberg et aLl l2008l : and NH3 at 9.0 /im and 



m 



Pontoppidan et al. 



2008 



CH,OH in BottineUi et al. 
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Fig. 2. — IRS spectra of the embedded sources in Serpens with rising mid-IR spectra. The 
prominent ice absorption features (at 6.0, 6.85, 7.7, 9.0 and 15.2 //m) are marked, together 
with the silicate absorption feature at 9.7 fj,m. The object # is indicated at the top right of 
each panel. 
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3.3. Disk Sources 



Young stellar objects with IR excess that have a flat or nega tiv e slope in the near - to 
mid-IR are identifled as disk sources. Following iGreene et al.lll994j and lEvans et al.ll2009l flat 
sources have —0.3 < a;2^m-24Atm < 0.3, Class II sources have —1.6 < Oi2^m-2i^lra < —0.3, and 
Class III sources have cx2^ra-2^^lr^^ < —1.6. Many of these sources also have silicate emission 
bands and a few show PAH features. In the following, these different features are analyzed 
and discussed in detail for the total of 94 disk sources. 



3.3.1. Disk Geometry: F30/-F13 and Cold Disks 

Disk geometry is inferred from the flux ratio between 30 and 13 yum (F30/-F13). High 
ratios (F30/-F13 > 15) yield edge-on disks, intermediate values (5 < F30/-F13 < 15 and 
1-5 ^ -FsoZ-^is ^ 5) identify cold disks and flared disks with considerable excess in the IR, 



Brown et al. 


2007; 


Merin et al. 


in prep.) 



). Figure [3] shows the distribution of F30/-F13 for the 



disk population in Serpens. 
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Fig. 3. — Distribution of the flux ratio between 30 and 13 /im, used as a proxy for disk 
geometry. 



Cold disks, sometime s referred to as transitional disks, as studied in lBrown et al.l (120071 ) 
and lMerfn et al.l ( lin prep. I ), present a peculiar spectral energy distribution (SED), with inner 
gaps or holes producing a region with no IR excess in the near- to mid-IR but a substan- 
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tial e xcess at mid- to far - IR wavelengths. Bin aries ( llreland &: Kraud 120081). planet forma- 



tion (IQuillen et al. 



2004 IVarniere et al.l l2006l ) and photoevaporation ( IClarke et al.l 12001 



Alexander et al.ll2006l ) are some of the suggested mechanisms for this removal of warm dust. 



Following iBrown et al.l ( 120071 ). -F30/-F13 is used to differentiate cold disks from the gen- 
eral disk population. Those wavelengths were chosen to avoid the silicate features and to 
probe the spectra steeply rising at mid- to far-IR following a deficit at shorter wavelengths, 
characteristic of cold disks. Eight objects (# 9, 21, 69, 82, 90, 113, 114 and 122) in this 
sample are classified as cold disks based on 5 < F30/F13 < 15, corresponding to 8.5% of 
the disk population. These spectra are shown in Figure HI For clarity, the original spectra 
are shown in grey, whil e in black a binned vers ion is overplotted. The binning was done 
using a Savitzky-Golay (ISavitzky &: Golaylll964l ) filter of order and degree 3, for illustra- 
tive purposes only. This is the most unbiased survey of cold disks to date, based entirely on 
spectroscopy in the critical mid-IR range, and should thus give a representativ e fraction of 



such disk s. These objects are also part of a larger cold disk sample, studied by iMerm et al. 



(jin prep. I ), where the authors have selected cold disks from photometric criteria in all five 
c2d clouds. 
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Fig. 4. — IRS spectra of the cold disks in Serpens. Object # is indicated in each panel. In 
grey are the original spectra, while in black a binned version is overplotted. The dotted line 
represents the stellar photosphere, for comparison. 



3.3.2. PAH Emission 

Polycyclic aromatic hydrocarbons (PAH) are large molecules whose emission is excited 
by UV radiation out to large distances from the star. Thus, these features are good diag- 
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nostics of the amount of stellar UV intercepted by the disk surface. 



Almost half of the s a mple of He rbig Ae/Be stars (young s tars o f intermediate masses) 
studied by iMeeus et al.l ( 120011 ) and lAcke &: van den Anckerl (120041 ) show PAH emission. 
Geers et al. ( 20061 ). on the other hand, show that PAH emission is not as common among 
low-mass YSOs as it is for more massive stars. The authors argue that the lack of features 
is not only due to their weaker UV radiation fields, but also to a lower PAH abundance 
compared with the general ISM. 
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Fig. 5. — IRS spectra of the disk sources which show clear PAH features. The dashed lines 
indicate the positions of the most promin ent bands (at 6.2, 7. 7, 8.6, 11.2, 12.8 and 16.4 /im). 
The object # and its spectral type (from lOliveira et al.ll2009l ) are indicated at the top right 
of each panel. Note that the ratios of different line strengths differ between sources. 
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In the current sample, 5 sources have clear PAH emission, presented in Figure [51 Al- 
though a great variation in shapes and strengths is seen, all sources show the features at 6.2, 
7.7, 8.6, 11.2 and 12.8 fim. Objects #52, 97 and 131 also show a feature that could be the 
PAH band at 16.4 fim, although the feature does not match perfectly in position. 



The central stars in these 5 sources were studied with optical spectroscopy in lOliveira et al 



(120091 ). Spectral types and luminosities were derived, yielding ages and masses for each ob- 
ject. Objects #52, 98, 120 and 131 have masses between 1.9 and 2.7 Mq (spectral type G3 
to A2), while object #97 has lower mass, of around 0.45 Mq (M2). Objects #97 and 98 are 
spatially very close to each other and have very similar PAH features. IRAC Band 4 (8 fim) 
images show that there is extended nebulosity surrounding both positions so that the PAHs 
are likely associated with this nebula rather than the disks. 

Discarding these 2 sources, this sample indicates that PAHs are present in 3.2% of the 



disk p opulation studied, which is below the 11-14% of PAH detection rate found by iGeers 
( 120071 ). This discrepancy could be due to the two samples being differently selected, with 
the current sample having a larger fraction of late-type stars. The percentage of PAHs 
present in disks in Serpens is, however, very comparable with that in Taurus (4% of the disk 
population, § 14. 3p . 

Interestingly, none o f the cold disks show PA Hs, different from expected if grain settling 
enhances PAH features (IDullernond et al.ll2007l) and observed for other cold disk samples 
studied which do show PAHs (IBrown et al.l 120071 : iMerm et al.llin prep. I ). In addition, the 



objects with PAH emission have standard -F30/-F13 flux ratios, falling in the 'flat' and 'flared' 
regions of Figure |3l 



Only object #120 shows a silicate feature (at 20 /zm, as discussed in § l3.3.3p . in addition 
to PAHs. This alludes to an anti-correlation between PAHs and silicates that could be simply 
explained in terms of contr ast, i.e., it is hard to detect PAH features when a strong 10 /im 
silicate feature is present (iGeers et al.ll2006l ). In this sense, PAHs should be more easily 
detected in disks without silicates, since the falling spectrum provides an increase in the 
feature to continuum ratio improving contrast. 



3.3.3. Silicate Emission 



S ilicate emission has been observed around numerous circumstellar disk s fe.g..lvan Boekel et al. 



2003; 


Przveodda et al. 


2003 


2008; 


Olofsson et al. 


2009). 



pared to primitive interstellar material. The strength and the shape of silicate features (at 
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10 and 20 fim) contain important information on the dust sizes in the surface layer of the 
disk. The 89 disk sources in this sample were inspected for their silicate characteristics. 
Their spectra can be seen in Figures [HI to [S] (completed with Figures S] and ED , arranged in 
order of decreasing strength of the 10 fim feature and decreasing a2-24^(m spectral slope. No 
extinction correction has been applied to these spectra (see § 13.3.51) . 




5 10 15 20 S5 30 35 5 10 15 SO S5 30 35 5 10 15 20 2b 30 35 



X (/Am) 



Fig. 6. — IRS spectra of the disk sources in Serpens. Object # is indicated in each panel. 
In grey are the original spectra, while in black a binned version (using a Savitzky-Golay 
(ISavitzky fc Golaylll964j ) filter of order and degree 3) is overplotted. 



Five of these sources (5% of the disk population) have spectra that are featureless, i.e.. 
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Fig. 7. — IRS spectra of the disk sources in Serpens, contd. 



they have neither the 10 nor the 20 fim sihcate features in emission, although they present IR 
excess (objects # 38, 42, 59, 83, and 87). It is found that when the 10 fim feature appears, 
the 20 yum band is also present (by visual inspection of the continuum-subtracted spectra, 
see §E33]). The contrary is not true, as objects #13, 20, 24, 32, 62, 64, 69, 70, 109, 116, 
127, 139, 14 3 and 145 have an ob s erved 20 fim feature but no apparent 10 /im feature (also 
reported bv iKessler-Silacci et all fl2006h for T Cha, SR21 and HD 135344B). Also, object 
#120 shows the 20 /im band and no 10 /im feature, as this region is dominated by PAH 
features in this source (see § I3.3.2p . When present, the 10 /im feature (S^l^^) is usually 
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Fig. 8. — IRS spectra of the disk sources in Serpens, contd. 



stronger than the 20 fim ("Spg^™) one as is common in disk sources. 



CrystaUine sihcate featu r es (as discussed by e.g. iBouwman et aLl2001 



2004 : Ivan Boekel et all 120041 : lApai et al.ll2005l : lOlofsson et al 



2009 



2008 



Sargent et ah 



Acke fc van den Ancker 
2OO9I ) are 



seen in several sources. A detailed analysis in terms of statistics and crystalline fraction of 
these features is non-trivial and beyond the scope of this paper, but will be subject of a 
future publication. 
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3.3.4- The Silicate Strength- Shape Relation 



As first shown for Herbig Ae stars (Ivan Boekel et al.ll2003l ). tliere is a dependence be- 
tween grain size and silicate sliape and strengtli for botfi 10 and 20 /xm features: as the sizes 
of the grains grow, the silicate features appear more and more flattened and shift to longer 
wavelengths. 

A statistical analysis of these two features can be done using the following quantities: 
^pcak is ^ measure of the strength of the feature centered at 10 /zm and S'n.a/S'g.s is a proxy 
for the shape of the feature (the smaller this value, the more peaked the feature is); S'^^^^ 
measures the strength of the 20 pm fe ature, while S'23.8/5'19 is a proxy for this feature shape. 
Following iKessler-Silacci et al.l (|2006[ ). Sx is defined as: 



1 

26 



X+5 



X-5 



d\. 



In this paper, we use 5 = 0.1 /im. 



Kessler-Silacci et al.l (|2006[ ) found a global trend (best fit to observations) for both 



features (Figure 9 in their paper and [9] in this one), indicating that the feature shape and 
strength are physically related, as expected if the silicate features are indeed tracers of grain 
size. The relationship is much tighter for the 10 ixm feature, as also seen by the great spread 
in the correlation for the 20 //m feature. 

Following their procedure, we first subtracted the continuum of each spectrum by fitting 
a second order polynomial to the following regions: blue-ward of the 10 /xm feature, between 
the 10 and the 20 /xm features, and red-ward of the 20 /xm feature (6.8-7.5, 12.5-13.5 and 
30-36 yum). The peak fluxes at 9.8, 10, 11.3, 19, 20 and 23.8 /xm were then determined 
(see Table |2]). Figure [9] shows t he results for our sample of disks, excluding those with PAH 
emission. The best fit found by iKessler-Silacci et al.l (120061 ) for the c2d IRS first look sample 
of YSOs has been i ndicated for compari s on (so lid line). The agreement of the Serpens data 
with the best fit of iKessler-Silacci et al.l ( 120061 ) is excellent, and more evident for the 10 /im 
than for the 20 /im feature due to the larger scatter in the latter group. 



Qlofsson et al.l ( 120091 ) generated synthetic 10 /xm features calculated for different grain 



sizes and compositions. Their models are generated for amorphous silicates of olivine and 
pyroxene stoichiometry and a 50:50 mixture, with grain size varying between 0.1 and ~6 /xm. 
Those models are overplotted in the left panel of Figure IHl with open symbols corresponding 
to different grain sizes. The comparison of the data presented here with these models shows 
that the majority of our sample lies in a region consistent with an opacity dominated by 




Fig. 9. — In the left panel, black dots show the ratio of peaks at 11.3 to 9.8 fim (5*11. 3/5*9. g) 
plotted against the peak at 10 fim (5'p|!^™). In the right panel, the ratio of peaks at 23.8 
to 19 /im (>S'^.-i.s/>S'ig) is plo tted against the peak at 20 /im {S^^i^^). The best fit found by 
Kessler-Silacci et al.l (120061 ) for the c2d IRS first look sample of YSOs (grey dots) has been 
indicated for comparison (solid black line). Green squares are the cold disks in this sample 
(see ^ 13.3.111. Colored c urves are derived from theoretical opacities for different mixtures 



by lOlofsson et al.l (|2009[ ). The open circles correspond to different grain sizes, from left to 
right 6.25, 5.2, 4.3, 3.25, 2.7, 2.0, 1.5, 1.25, 1.0 and 0.1 fim. Typical uncertainties for peak 
strength are ~0.1, while typical uncertainties for the peak ratios are ~0.12. 



grains with sizes larger than 2.0 /im, having no features consistent with grain sizes smaller 
than 1.5 fim. A considerable fraction of the sample is consistent with grains as large as 6 
/im (the largest grain size modeled). The precise sizes depend on composition and treatment 
of the opacities, but sizes of a few /im are clearly implicated. It is important t o note that 
crysta llinity also affects th e shap e of silicate features. However, as shown by lApai et al. 
(I2OO5I ) and lOlofsson et al.l ( 120091 ) . the effect of crystallinity is orthogonal to that of grain 
sizes in the strength versus shape plot (Figure [9]). Grain sizes are found to be the dominant 
parameter changing the shape of the 10 /x m silicate feature, whereas crystallinity introduces 
scatter (Figure 13 of lOlofsson et al.ll2009l ). 



All cold disks present silicate features in emission. According to the models of grain 
sizes (left panel of Figure |9]), the cold disks in this sample have grains bigger than 2.7 fim, 
with the bulk of the sample presenting grains as big as 6.0 /im, according the these models. 
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3.3.5. Effects of Extinction 



The amount of foreground extinction is, for some sources, substantial enough to affect 
the appearance of sihcate emission features in the mid-infrared. Because the extinction law 
has strong resonances from silicates, the strength and shape of any silicate emission feature 
may be significantly affected if a large column of cloud material is present in front of a given 
source. Here, the potential effects of extinction on statistical results, such as those presented 
in Figure [H are discussed. 



The model presented by IWeingartner fc Draind (120011 ) for an absolute to selective ex- 
tinction of Ry = 5.5 has been found to be a reasonable match to the observed dark cloud 



extinction law f IChapman et al.l 120091 : iMcClurd l2009l : IChapman fc Mundvl 120091 ) . and it is 
assumed that this law holds for Serpens. One caveat is the lack of resonances due to ices 
in this model, but given other uncertainties in the determination of extinction laws, this is 
probably a minor contribution for the purposes of this work. 



The opacity caused by a silicate resonance is defined as (e.g. lDraindl2003l ): 



Ata = Ta - Tc, 



(2) 



where tq is the optical depth of the "continuum" opacity without the silicate resonance. 
The relation between optical extinction and Ar is, using the fact that the optical depth is 
proportional to the extinction coefficient C*^^*: 



A, 



/^cxt 



Ata 0.921ArA 



0.921ACft 



= Xx mag. 



(3) 



Using Equations [T] and |3] and by approximating 5 — )■ 0, the silicate strength corrected 
for extinction is: 



Sx 



exp(-ArA) exp{-Av/Xx 



(4) 



For the Ry = 5.5 extinction law of IWeingartner fc Draind (120011 ). the relevant values for Xa 
are 19.6, 20.9 and 41.6 mag for 9.8, 10.0 and 11.3 fim, respectively. 

Using this relationshi p, it is possible to c orrect the position of sources with known 
values for Ay in Figure [9l lOliveira et al.l (120091) measured Ay values for 49 of the 89 disks 
studied here albeit with some uncertainty. The resulting extinction-corrected strength-shape 
distribution is shown in Figure [10] (in green) and compared with the uncorrected distribution 
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(in black). It is seen that extinction moves points along curves that are almost parallel to 
the relation between shape and strength. Most points, having Ay < 5 mag, are not changed 
enough to alter the slope of the relation, and the median strength is unaffected. However, a 
few highly extincted disks are corrected by large amounts, as can be seen by the extinction 
arrows in Figure [TOl 




1.0 1.5 2.0 2.5 3.0 3.5 

peak 



Fig. 10. — Comparison between the observed values (in black) and the extinction corrected 
values (in green) for the 10 fim silicate feature. The arrows give an indication of the effects 
of extinction on the feature, however their length depends on the starting point (see text for 
details) . 



There are 40 YSOs in the sample for which lOliveira et al.l (120091 ) do not estimate Ay- 
Some of these may be highly extincted {Ay ^ 10 mag) sources (not being bright enough for 
optical spectroscopy using a 4-m class telescope). The correction for extinction, as derived 
above, would be greater for such objects, possibly introducing a displacement between the 
distributions of extinction corrected and uncorrected features. However, it is important to 
note that this assumes that the Ry = 5.5 extinction law is valid for the densest regions of 



molecular clouds. IChiar et al.l (120071 ) recently found that at Ay > 10 mag, the relationship 
between optical depth and extinction derived for diffuse ISM is no longer valid, possibly 
due to grain growth, and t hat X Q 7^,,rn is significantly larger in this regime. If true, then 
the IWeingartner fc Draind (120011 ) extinction law represents a "worst case" scenario. Since 
most of the measured extinctions are small, and given the uncertainty in extinction laws and 
the lack of Ay data for a significant fraction of the sample, the uncorrected silicate feature 
strengths are used in the remainder of the paper. 
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3.3.6. Statistical analysis of silicate features 



Figure [TT] shows the strengths of the observed 10 and 20 /xm features with no apparent 
correlation. To quantify a correlation, a Kendall r rank correlation coefficient is used to 
measure the degree of correspondence between two populations and to assess the significance 
of said correspondence. In other words, if there is a correlation (anti-correlation) between 
two datasets, the Kendall r rank coefficient is equal to 1 (-1). If the datasets are completely 
independent, the coefficient has value 0. For the strengths of the 10 and 20 /im silicate 
features, r = 0.26, meaning at best a weak correlation. 

The lack of apparent correlation between the 10 and the 20 /im features, in both strength 
and appearance, suggests that they are emitted by different regions in the disk. Indeed, sev- 
eral authors have shown that the 10 /im feature is emitted by warm dust in the inner region, 
while the 20 Atm feature originates from a colder component, further out and deeper into the 



disk (e.g. , iKessler-Silacci et al.ll2007l : iMerm et al.ll2007l : iBouwman et al.ll2008t lOlofsson et al. 



20091 ). Therefore, the absence of the 10 ixm feature but presence of the 20 /xm band for a 
given source implies that such disk lacks warmer small dust grains close to the star. This 
reminisce the disk around th e Herbig Ae star HD 100453 studied by lMeeus et al.l (120021 ) and 
Vandenbussche et al.l ( 120041 ). who find that the absence of the 10 /xm silicate feature can be 
fitted with a model deprived in warm small grains. 16% of the disk sources shown in this 
sample fit this scenario. A possible explanation for such a scenario is that those disks have 
holes or gaps on such a small scale (< 1 AU) that the holes do not produce a strong signature 
in the spectra probed by our data, as the cold disks do. 



4. Discussion 



In this section we discuss the properties of the YSO sample presented in sections 
13.21 and 13.31 with respect to environ ment, and compar e the results with those for another 
nearby star- forming regioii , Taurus (lFurla.n et al.ll2006[ ). as well as the full c2d IRS sample 
dKessler-Silacci et al.lbooel : loiofsson et al.lboOQl ). 



4.1. Cluster vs. Field Population 

Comparison of the disk properties between the cluster and field populations determines 
the importance of environment in the evolution of these systems. If the evolution of the 
disks in clusters is found to follow a different pace than that for more isolated stars, this can 
set important constraints on disk evolution theories. 
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Fig. 11. — Relative strength of the 10 vs. 20 yum sihcate features. Typical error bars are 
shown in the bottom right corner. No obvious correlation is found. 



As discussed in lHarvey et al.l (j2006h . the brightest YSOs in Serpens appear to be concen- 
trated in clusters, but a more extended young stellar population exists outside these clusters. 
For the determination of clusters and their boundaries, we follow the method developed by 
J0rgensen et al.l (120081 ) for Ophiuchus and Perseus. Using volume density and number of 
YSOs as criteria, the method consist s of a nearest neighb or algorithm, and has as input the 
complete sample of YSOs in Serpens ( iHarvey et al.ll2007bl ). The associations are divided into 
loose (volume density of 1 Mopc"^, blue contour in Figure fT2l) and tight (volume density of 
25 MqPc~^, yellow contours in Figure [T^ . Furthermore, the associations are divided into 
clusters (more than 35 members) and groups (less than 35 members). The results of this 
method in Serpens yields 2 clusters (A and B) and 3 groups (C, D and E), that can be seen 
in Figure [121 Individual memberships are marked in Table El With the exception of one, 
virtually all YSOs in Serpens are within the blue contour, i.e., all YSOs be long to, at least 
a loo se association. This differs, for instance, from Ophiuchus and Perseus (jj0rgensen et al. 



20081 ). where s ome YSOs are cona pletely isolated. For number statistics of clustering in all 
c2d clouds, see lEvans et al.l ( 120091 ). 



Table El lists the number of objects in each tight association. It appears that the ratio 
of disk to embedded sources in a given cluster increases with distance from the densest part 
of the cluster, where Group C is located (see Figure [T2|) . 
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Fig. 12.— C 



J0rgensen et al. 



usters and groups in Serpens as defined according to the criteria in 
(120081) . The volume density contours are overlaid on the Se rpens extinc- 



tion m ap (lEvans et al.l 120071 ). The red dots are the YSOs in Serpens from iHarvey et al. 



(l2007bl ). Black contours indicate volume densities of 0.125, 0.25, 0.5, 2.0 and 4.0 MqPc"^, 
blue contour corresponds to volume density of 1 MqPc"^ and yellow contours to a volume 
density of 25 MqPc^^. 



To compare the populations, all objects belonging to any of the tight associations were 
grouped into the "cluster population", for better number statistics. Objects not belonging 
to any cluster were grouped into the "field population" . Figure [13] shows the comparison 
between these two populations for 3 quantities derived from the IRS spectra. In the left 
panel, the flux ratio between 30 and 13 fim (F30/-F13) is an indication of disk geometry. The 
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middle and right panels show the peak intensity of the silicate features at 10 and 20 /im, 
respectively. For all 3 quantities, the two populations (cluster and field) are statistically 
indistinguishable. A two sample Kolmogorov-Smirnov test (KS-test) was performed for each 
quantity and the results show that the null hypothesis that the two distributions come from 
the same parent population cannot be rejected to any significance (18%, 13% and 79% for 
i^3o/i^i3, ^peaiT ^nd ^p"^^, respectively). 
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Fig. 13. — Comparison between the clustered (dot-dashed black line) and field popula- 
tions (solid blue lin e) of Serpens, with data from the c2d IRS program (dashed red line. 



Olofsson et al.ll2009l ). In the left panel, the flux ratio between 30 and 13 /im (-F30/-F13) is 
an indication of the disk geometry. The middle and right panels show the strength of the 
silicate features at 10 and 20 fim, respectively. 

Differences between the cluster vs. field populations will be further investigated with 
ancillary data (e.g., relative stellar ages and masses) and modeling (e.g., disk sizes). How- 
ever, the IRS spectra allow the conclusion that no significant differences are found for disk 
geometry or the grain size distribution in the upper layers of circumstellar disks in clustered 
compared with field stars. The latter result only applies to the inner disk, as traced by 
silicate features. The outer disk may still be different between cluster and field populations. 



4.2. Comparison with Other Samples 



The young stars with disks observed by the c2d team with IRS s pectroscopy (c2d IR S 
sample) are scattered across the sky in the 5 molecular clouds studied (lOlofsson et al.ll2009l ). 
The other 4 clouds studied by the c2d are Chamaeleon II, Lupus, Perseus and Ophiuchus. 
All clouds are nearby (within 300 pc), span a range of star-formation activity, have typical 
median ages of a few Myr and have a spread betwee n more (Perseus a nd Ophiuchus) and 
less (Cha II and Lupus) clustered YSO populations (lEvans et al.ll2009l ). The results from 
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this sample are also compared to the results in Se rpens in Figure s [9] an d [T3l A conspicuous 
similarity is seen between the samples. Note that lOlofsson et al.l (120091 ) do not analyze the 
20 fjm silicate feature in this same manner and, therefore, S'^^^^ for the c2d IRS sample is 
missing from the right panel in Figure [T3i 



o 




3.5 



peak 



Fig. 14. — The flux ratio between 30 and 13 fim (F30/-F13) plotted against the peak at 
10 Atm (S^^'^^, black dots), compared against the T Tauri FEPS sample (red triangles, 
Bouwman et al.ll2008l ) and the c2d IRS sample (blue stars, lOlofsson et al.ll2009[ ). Typical 
uncertainties for the 10 fim peak strength are ~0.1, while typical uncertainties for the flux 
ratio are ~0.12. 



The models of iDuUemond fc DominikI ( l2008l ) conclude that, if sedimentation is the 
unique reason for the variety of observed strength and shape of the 10 /im silicate fea- 
ture, then this feature is strong for weak excess in mid- to far-IR (as probed by F30/F13), 
and vice-versa. However, when studying a small sample of T Tauri stars from the Spitzer 
Legacy Program "The Formation an d Evolution of Planeta ry Systems: Placing Our Solar 
System in Context" (FEPS) sample, iBouwman et al.l (|2008[ ) found the opposite: a trend in 
which weak F30/F13 correlates with a weak feature. A confirmation of this trend for a larger 
sample implies t hat sedimentation alone c annot be the sole cause for the diversity of observed 
silicate features. iDuUemond fc DominikI (120081 ) argue that dust coagulation must then play 
a vital role in producing different silicate profiles. In Figure HU we populate this diagram 
with the Serpens sample (black dots) and the c2d IRS sample (blue stars), as well as the 
Bouwman FEPS sample (red triangles). This large combined sample shows no correlation 
between the strength of the 10 /im silicate feature and F30/-F13 (r = 0.07) and therefore does 
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not support eit her the correlation (as s e en bylBouwman et al.ll2008l ) or the ant i- correlation 
(as modeled by iDuUemond &: DominikI ( 120081 ) for sedimentation alone) between the IR flux 
excess and the strength of the 10 /im silicate feature. 



4.3. Comparison with Taurus 



The Taurus Molecular Cloud is the best characterized star-forming region to date, due 
to its proximity and relatively low extin ction. With young stars arid their surrounding disks 



studied for more than two decades (e.g.. iKenyon fc Hartmannlll987l . ll995[ ). its members have 



been well characterized at a wide range of wavelengths. Taurus has thus become the refer- 
ence for comparison of star-forming regio ns. Compared to Serp ens, Taurus seems to be a 



some what younger (2.0 Myr median age, iHartmann et al.l 12001 



vs. 4.7 Myr lOliveira et al. 



20091 ) and has low er star forming rate . Due to uncertainties in pre-ra ain sequence age de 



terminations (e.g., iBaraffe et al.ll2009l : iHillenbrandl l2009l : iNaylorl |2009| ). this difference may 



not be signiflcant. However, cluster ages are statistically more likely to be different than the 
same. 

In a campaign similar to that presented here, IRS spectra were obtained for a sample 
of 139 YSOs in Taurus, as part of a larger I RS guaranteed - time observing program. These 
data were presented in t wo separate papers: iFurlan et al.l ( 120061 ) treated the disk sources, 
while iFurlan et al.l ( l2008l ) presented the embedded population in Taurus. 



Because the young stellar population in Serpens has been discovered and characterized 
using IR observations, this sample does not include young stars withou t disks (and therefore 
no IR excess, also called Class III sources, as deflned in iLadal ( 1l987l )). For this reason, a 
comparison between the Class III sources of the two clouds does not make sense. 26 Class III 
objects were studied in Taurus. These objects present featureless spectra with very little IR 
excess at longer wavelengths. It is, however, possible to compare the IR excess populations 
of both regions. Out of the entire YSO sample studied in Taurus, embedded sources (or 
Class I) amount to 20% while in Serpens they amount to 18%, a very comparable number. 
This also matches the typical percentage of cl ass I sources in the global YSO statistics of the 
flve clouds from the c2d photometric sample (lEvans et al.ll2009l ). In Serpens only 5 objects 
have featureless spectra. 

Excluding the embedded sources to analyze the disk population alone, objects having 
both 10 and 20 /xm silicate features amount to 72% in Taurus, comparable to 73% of the 
Serpens disk population. Also comparable is the percentage of disks with PAH emission: 3% 
in Serpens and 4% in Taurus. 
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In both regions, all disks with a 10 /im silicate feature also show the 20 fim one. However, 
a difference is found for disks with only the 20 /im feature in emission: they amount to 17% 
of the disk population in Serpens and only 4.5% in Taurus. These statistics are summarized 
in Table H 
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Fig. 15. — Top: Distribution of the 10 /im peak strength for Serpens (black) and Taurus 
(red). Bottom: Strength and shape of the 10 fim silicate feature for both Serpens (black) 
and Taurus (red). 



IRS data on the 85 disks sources presented in iFurlan et al.l (120061 ) were obtained from 
the Spitzer archive and reduced in the same manner as our Serpens data. The same method 
described in Section 13.3.41 was applied to these data in order to compare the processes 
affecting the dust in both regions. Figure [15] (bottom) shows the strength versus shape of 
the 10 fim silicate feature for both Serpens (black dots) and Taurus (red dots), as well as the 
distribution of peak strength for both regions (top). Serpens and Taurus present remarkably 
similar distributions, with populations clustered around flatter features and bigger grains, 
and almost no strongly peaked silicate emission sources. A KS-test shows that the hypothesis 
that the distributions are drawn from the same population cannot be rejected (82%). 
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4.4. Implications 

Comparisons between disks around stars tliat were formed in clusters or in isolation 
in the Serpens Cloud, explored in § 14.11 indicate very similar populations. Thus, statisti- 
cally no differences are found in terms of both disk geometry (probed by the ratio F30/-F13) 
and processes affecting the dust (probed by the silicate features). This suggests that local 
environment does not affect the evolution path and timescale of disks. 

Even more remarkable is the agreement between the YSO populations of Serpens and 
Taurus, as well as the c2d IRS sample spread over 5 clouds, in terms of silicate features, 
as shown in § 14.21 and 14.31 Even though there are clear differences in silicate features from 
source to source within a cloud, the overall distribution of feature shapes is statistically 
indistinguishable: in all three samples, each containing of the order of 100 disks, the bulk of 
the sources has a rather flat silicate profile, and a tail of peaked shapes. 

If the difference in median ages is significant, the similarity seen in Figure [TS] indicates 
that a 2-3 Myr difference in age is not reflected in a concurrent evolution of the average disk 
surface dust properties. This indicates that the dust population in the disk surface is an 
equilibrium between dust growth and destruction processes, which is maintained at least as 
long as the disk is optically thick. 

The process of grain growth through coagulation and settling to the midplane has been 



shown to be much too short to be consistent with disk observations (jWeidenschillingl Il980 



DuUemond &: Dominikll2005l ). That means that the small grain population must be replen- 
ished somehow, which could happen by fragmentation of bigger aggregates and turbulent 
mixing. It is widely accepted in the literature that the 10 /xm silicate fea ture is representative 



of th e dust in the surface layers of the disk at a few AU from the star (IKessler-Silacci et al. 



20071 ). Significant evolution could still take place in the disk mid-plane, which is not traced 



with these observations. Indeed, millimeter observations have indicated that disk mid- 



planes are abundant in grains with mm/cm sizes (INatta et al.ll2004l : iRodmann et al.l 12006 



Lommen et al.ll2007l ). That means that, if this population of small grains in the disk surface 
must be replenished by fragmentation of bigger grains followed by turbulent mixing, the 
surface dust is an indirect tracer of midplane grains. 

If the age difference between Serpens and Taurus is significant, then there exists a sta- 
tistical equilibrium of processes of dust coagulation and fragmentation in the disk lasting 
at least a few million years that is independent of which YSO population is being stud- 
ied. The observable IR dust characteristics of a cluster population of protoplanetary disks 
does not appear to depend on cluster properties (age, density) within the first 5 Myr, for 
relatively low-mass clusters. No specific property or event that may influence the infrared 
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dust signatures of a single disk (stellar luminosity, spectral type, presence of a companion, 
disk-planet interactions, disk instabilities etc.) produces detectable temporal or spatial evo- 
lutionary effects visible in the distribution of disk properties of an entire cluster within this 
time frame. 

This is consistent with a picture in which a specific disk may change its appearance on 
short time scales (much less than 1 Myr), but in a reversible way. That is, the properties o f 



disk surfaces may oscillate between different states ( iBary et al.ll2007l : iMuzeroUe et al.ll2009l ). 
If some evolution of the disk surface is irreversible, and happens at a given rate, a stable 
distribution of surface properties would not be seen over time. This requires that the effect 
on the disk surface of any reprocessing events has to be erased. 

Another option is that the disk surface properties are determined by a single parameter, 
such as the initial conditions of the formation of such disks. In this scenario those properties 
should be kept "frozen" over the observed time scale of ~few Myr. However, this cannot be 
the case for objects like EX Lup, where real time changes of the disk surface properties after 
the initial formation of the disk have been observed on time scales of just a fe w years (namely 



cryst allization through thermal annealing of the dust in the disk surface, [Abraham et al. 



20091 ). 



For this theory to work, the putative oscillation of states proposed here must be stable 
over disk lifetimes of ~5 Myr. As long as the disks are gas-rich and optically thick, as in the 
three samples studied here, we do not see a distinguishable evolution. It will be important 
to search for surface evolution indicators in even older clusters. Such comparisons will yield 
constraints on the global time scale of disk evolution, which in turn will help constrain the 
importance of the processes that play a vital role in the dissipation of the disks. 

The scenario in which a considerable population of small particles is still present after 
bigger aggregates are formed and altered through multiple events is consistent with evidence 
from primitive chondrites in our own Solar System. Chondritic meteorites are observed to 
contain fine-grained dust-like matrices formed after coarse-grained materials such as chon- 
drules and calcium-aluminum-rich inclusions (CAIs). This indicates that the solar nebula 
underwent violent reprocessing events in the feeding zones of parent bodies at 3-4 AU (see 



Pontoppidan fc Brearleyi ( l2010l ) for a recent review on Solar System dust in an astrophysical 
context). These events included systemic melting of large dust aggregates and possibly even 
evaporation and re-condensation of silicate grains. In fact, presolar material is present in 



prim itive chondrites only at the trace level: ~100 ppm (ILodders fc Amarill2005l : IZinner et al. 



20071 ). testifying to the complete evaporative destruction of ISM dust at a few AU in proto- 



planetary dust. 
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5. Conclusions 

We present Spitzer/IRS spectra from a complete and unbiased flux limited sample of 
IR excess sources found in the Serpens Molecular Cloud, following the c2d mapping of this 
region. 

• Among our total of 147 IRS spectra, 22% are found to be background contamination 
(including background stars, redshifted galaxies and a planetary nebula candidate). 
This high number is not surprising given the position of Serpens, close to the galactic 
plane. 

• Excluding the background objects from the sample, the bona flde set of YSOs amounts 
to 115 objects. The embed ded to disk source ratio is 18%, in agreement with the ratios 



derived from photometry (lEvans et al.l 120091 ) for the flve c2d clouds. 



Disks with PAH in emission amount to 2% of the YSO population, but 3% of the disk 
population. Only G and A star show PAH emission. In the disk population, 73% 
show both silicate emission features, at 10 and 20 /im, while 17% show only the 20 
fim feature. 4% of the disks show featureless mid-IR spectra. Only one source, #120, 
shows both PAH and silicate in emission. 

Our YSO population in Serpens is very similar to that in Taurus, also studied with 
IRS spectra. In both regions about 70% of the disk sources present both 10 and 20 
/xm silicate features in emission and the 10 fim feature is never seen without the 20 fim 
feature. The only significant difference between the two populations is in the number 
of sources lacking the 10 /im feature but showing the 20 fim feature. This may be 
indicative of small holes (< 1 AU) in these sources. 

The silicate features in the IRS spectra measure the grain sizes that dominate the mid- 
IR opacity. The relationships between shape and strength of these features present 
distributions very similar to th ose obtained from ot her large samples of young stars. 



Comparison with the models of lOlofsson et al.l ( 120091 ) yield grains consistent with sizes 
larger than a few fim. 

No significant differences are found in the disk geometry or the grain size distribution 
in the upper layers of circumstellar disks (probed by the silicate features) in clustered 
or field stars in the cloud, indicating that the local environment where a star is born 
does not influence the evolution of its harboring disk. 

Quantitatively, the shape and strength of the 10 /im silicate feature were used to 
compare both Serpens and Taurus as well as a large sample of disks across flve clouds. 
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indicating remarkably similar populations. This implies that the dust population in the 
disk surface results from an equilibrium between dust growth and destruction processes 
that are maintained over a few million years. 
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Table 1. Observed objects in Serpens 





ID'' 


c2d ID (SSTc2dJ) 


Q;2(im — 24^m'^ 


Program'' 


AOR 


Obs. Date 




Classification 


1 


1 


18275383-0002335 


-1.42 


G03 


17891072 


2007-04-19 


1.31 


Disk 


2 


2 


18280503+0006591 


-2.29 


G03 


17883648 


2007-04-24 


0.35 


BG star 


3 


3 


18280845-0001064 


-0.70 


G03 


17883136 


2007-04-24 


1.52 


Disk 


4 


4 


18281100-0001395 


-2.13 


G03 


17888256 


2007-04-29 


0.27 


BG star 


5 




18281315-1-0003128 


-1.84 


c2d 


13210368 


2005-05-20 


0.27 


BG star 


6 


5 


18281350-0002491 


-1.08 


G03 


17883136 


2007-04-24 


1.71 


Disk 


7 


6 


18281501-0002588 


-0.05 


G03 


17883136 


2007-04-24 


2.54 


Disk 


8 


7 


18281519-0001407 


-1.25 


G03 


17889536 


2007-04-30 


2.43 


Disk 


9 


8 


18281525-0002434 


-0.85 


G03 


17882368 


2007-04-25 


5.37 


Cold disk 


10 


9 


18281629-0003164 


-1.13 


G03 


17884160 


2007-04-19 


1.83 


Disk 


11 




18281757-0006474 


0.65 


G03 


17890560 


2007-04-28 


3.36 


Galaxy 


12 




18281757-1-0016065 


0.53 


G03 


17888768 


2007-05-05 


4.72 


Galaxy 


13 


11 


18281981-0001475 


-0.89 


G03 


17889792 


2007-04-30 


4.09 


Disk 


14 


13 


18282143+0010411 


-1.44 


G03 


17888000 


2007-04-30 


1.38 


Disk 


15 


14 


18282159+0000162 


-0.99 


G03 


17883136 


2007-04-24 


2.56 


Disk 


16 


15 


18282432+0034545 


-2.01 


G03 


17887744 


2007-04-24 


0.32 


BG star 


17 




18282720+0044450 


1.82 


G03 


17884928 


2007-04-29 


11.04 


PN candidate 


18 


16 


18282738-0011499 


-1.98 


G03 


17883136 


2007-04-24 


0.32 


BG star 


19 


17 


18282741+0000239 


-1.97 


G03 


17887744 


2007-04-24 


0.43 


BG star 


20 


18 


18282849+0026500 


-1.12 


G03 


17890048 


2007-04-29 


1.13 


Disk 


21 


19 


18282905+0027560 


-1.06 


G03 


17884928 


2007-04-29 


10.38 


Cold disk 


22 




18283000+0020147 


-1.60 


c2d 


13210624 


2005-05-20 


0.29 


BG star 


23 




18283736+0019276 


-2.19 


G03 


17887744 


2007-04-24 


0.27 


BG star 


24 


20 


18284025+0016173 


-0.96 


G03 


17890048 


2007-04-29 


1.20 


Disk 


25 


21 


18284053+0022144 


-1.91 


G03 


17887744 


2007-04-24 


0.42 


BG star 


26 


22 


18284187-0003215 


0.53 


G03 


17883136 


2007-04-24 


3.39 


Emb 


27 


23 


18284403+0053379 


0.70 


G03 


17882880 


2007-04-28 


3.37 


Emb 


28 


24 


18284479+0051257 


1.04 


G03 


17885184 


2007-04-24 


10.18 


Emb 


29 


25 


18284481+0048085 


-1.01 


G03 


17888512 


2007-05-05 


2.76 


Disk 


30 


27 


18284497+0045239 


-1.35 


G03 


17886720 


2007-04-24 


0.77 


Disk 


31 


38 


18284559-0007132 


-0.77 


G03 


17882368 


2007-04-25 


0.51 


Disk 


32 


29 


18284614+0003016 


-0.64 


G03 


17889024 


2007-04-25 


0.66 


Disk 


33 




18284632-0011103 


0.53 


G03 


17890560 


2007-04-28 


11.17 


Galaxy 


34 




18284828-0005300 


-2.43 


G03 


17888256 


2007-04-29 


0.47 


BG star 


35 




18284938-0006046 


-2.01 


c2d 


13210624 


2005-05-20 


0.28 


BG star 


36 


32 


18285020+0009497 


-0.24 


c2d 


13461505 


2006-04-31 


1.98 


Disk 


37 


33 


18285039-0012552 


-1.40 


G03 


17884160 


2007-04-19 


0.32 


BG star 


38 




18285060+0007540 


-2.31 


c2d 


13460736 


2005-09-09 




Disk 


39 


34 


18285122+0019271 


0.45 


G03 


17883904 


2007-04-29 


8.98 


Emb 


40 


36 


18285249+0020260 


-0.15 


G03 


17883136 


2007-04-24 


1.46 


Disk 


41 


37 


18285276+0028466 


0.07 


G03 


17887232 


2007-04-27 


9.50 


Disk 


42 


38 


18285362+0019302 


-0.64 


G03 


17887488 


2007-04-25 


2.09 


Disk 


43 


39 


18285395+0045530 


-1.16 


G03 


17886208 


2007-09-30 


1.96 


Disk 


44 


40 


18285404+0029299 


1.35 


c2d 


13461249 


2006-04-20 


18.21 


Emb 


45 


42 


18285450+0029520 


0.15 


c2d 


13460993 


2006-04-19 


35.49 


Emb 


47 


43 


18285489+0018326 


0.88 


G03 


17883392 


2007-04-29 


7.58 


Emb 
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48 


44 


18285529+0020522 


-0.29 


G03 


17884160 


2007-04-19 


1.43 


Disk 


49 


45 


18285580+0029444 


1.81 


G03 


17883392 


2007-04-29 


99.75 


Emb 


50 


46 


18285660+0030080 


1.84 


G03 


17889792 


2007-04-30 


21.37 


Emb 


51 


47 


18285719+0048359 


-0.32 


G03 


17888768 


2007-05-05 


11.27 


Emb 


52 


48 


18285808+0017244 


-2.11 


G03 


17888000 


2007-04-30 


1.23 


Disk 


53 


49 


18285860+0048594 


-1.06 


G03 


17887488 


2007-04^25 


1.55 


Disk 


54 


50 


18285946+0030029 


-0.59 


c2d 


13461249 


2006-04^20 


0.62 


Disk 


55 


51 


18290025+0016580 


-1.01 


G03 


17884416 


2007-04-30 


1.50 


Disk 


56 


52 


18290057+0045079 


-0.94 


G03 


17886208 


2007-09-30 


3.22 


Disk 


57 


53 


18290082+0027467 


-0.36 


G03 


17885440 


2007-04-24 


3.41 


Disk 


58 


54 


18290088+0029315 


-0.37 


c2d 


13210112 


2005-04-17 


1.82 


Disk 


59 


55 


18290107+0031451 


-0.49 


c2d 


13461249 


2006-04-20 


1.14 


Disk 


60 


56 


18290122+0029330 


-0.39 


c2d 


13461505 


2006-04^21 


1.58 


Disk 


61 


58 


18290175+0029465 


-0.84 


c2d 


13461249 


2006-04^20 


2.30 


Disk 


62 


59 


18290184+0029546 


-0.77 


c2d 


13210112 


2005-04-17 


1.32 


Disk 


63 


60 


18290210+0031210 


0.29 


G03 


17887232 


2007-04-27 


32.23 


Emb 


64 




18290215+0029400 


-0.86 


G03 


17885184 


2007-04-24 


2.89 


Disk 


65 


61 


18290286+0030089 


-0.23 


c2d 


13460993 


2006-04-21 


2.43 


Disk 


66 


62 


18290393+0020217 


-0.82 


G03 


17883136 


2007-04-24 


1.19 


Disk 


67 


63 


18290436+0033237 


0.95 


c2d 


13461505 


2006-06-21 


2.54 


Emb 


68 




18290442-0020018 


0.52 


G03 


17890560 


2007-04-28 


5.61 


Galaxy 


69 


64 


18290518+0038438 


-1.24 


G03 


17889792 


2007-04-30 


5.64 


Gold disk 


70 


65 


18290575+0022325 


-2.12 


G03 


17886976 


2007-04-29 


1.76 


Disk 


71 


66 


18290615+0019444 


-1.52 


G03 


17888512 


2007-05-05 


2.49 


Disk 


72 


67 


18290620+0030430 


1.69 


c2d 


13461249 


2006-04^21 


17.07 


Emb 


73 


68 


18290680+0030340 


1.67 


c2d 


13460993 


2006-04^21 


75.97 


Emb 


74 


69 


18290699+0038377 


-0.63 


G03 


17884416 


2007-04-30 


2.90 


Disk 


75 


71 


18290765+0052223 


-1.01 


G03 


17888512 


2007-05-05 


2.01 


Disk 


76 


72 


18290775+0054037 


-1.06 


G03 


17888000 


2007-04-30 


1.34 


Disk 


77 


73 


18290808-0007371 


-2.22 


G03 


17888256 


2007-04-29 


0.27 


BG star 


78 




18290843-0026187 


-2.28 


G03 


17885696 


2007-04-28 


0.31 


BG star 


79 


75 


18290910+0031300 


0.24 


G03 


17884928 


2007-04^29 


51.02 


Emb 


80 


76 


18290956+0037016 


-1.13 


G03 


17882368 


2007-04^25 


0.51 


Disk 


81 


77 


18290980+0034459 


-0.97 


c2d 


13210624 


2005-05-20 


0.85 


Disk 


82 


78 


18291148+0020387 


-1.34 


G03 


17887232 


2007-04-27 


11.90 


Cold disk 


83 


79 


18291249+0018152 


-1.06 


G03 


17889536 


2007-04-30 


2.06 


Disk 


84 


81 


18291407+0002589 


-2.19 


G03 


17885696 


2007-04-28 


0..30 


BG star 


85 




18291477-0004237 


-2.29 


c2d 


13210112 


2005-04-17 


0.24 


BG star 


86 


83 


18291508+0052124 


-1.89 


G03 


17882880 


2007-04^28 


0.45 


Disk 


87 


84 


18291513+0039378 


-1.50 


G03 


17889024 


2007-04^25 


1.01 


Disk 


88 


85 


18291539-0012519 


-1.72 


G03 


17885952 


2006-10-23 


0.50 


Disk 


89 


86 


18291557+0039119 


-0.56 


G03 


17884416 


2007-04-30 


2.51 


Disk 


90 


87 


18291563+0039228 


-0.89 


G03 


17885184 


2007-04-24 


6.52 


Cold disk 


91 


88 


18291617+0018227 


0.45 


c2d 


13210112 


2005-04-17 


2.11 


Emb 


92 


89 


18291969+0018031 


-1.23 


G03 


17885440 


2007-04-24 


2.10 


Disk 


93 


90 


18292001+0024497 


-2.04 


G03 


17886720 


2007-04-24 


0.29 


BG star 
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94 




18292050+0047080 


1.29 


G03 


17888768 


2007-05-05 


25.76 


Disk 


95 


93 


18292094+0030345 


-1.73 


G03 


17882880 


2007-04-28 


0.40 


BG star 


96 


94 


18292184+0019386 


-0.99 


G03 


17886976 


2007-04-29 


1.54 


Disk 


97 




18292250+0034118 


-1.65 


G03 


17889280 


2007-04-24 


0.29 


Disk 


98 




18292253+0034176 


-2.48 


G03 


17889024 


2007-0+25 


1.34 


Disk 


99 


95 


18292616+0020518 


-2.26 


G03 


17883648 


2007-04^24 


0.37 


BG star 


100 


96 


18292640+0030042 


-0.98 


G03 


17888768 


2007-05-05 


1.89 


Disk 


101 


97 


18292679+0039497 


-0.60 


G03 


17885440 


2007-04-24 


4.13 


Disk 


102 


99 


18292736+0038496 


0.29 


G03 


17886464 


2007-04-29 


3.95 


Emb 


103 


100 


18292824-0022574 


-0.41 


c2d 


13210368 


2005-05-20 


1.17 


Disk 


104 




18292833+0049569 


-0.48 


G03 


17887488 


2007-04-25 


1.85 


Disk 


105 




18292864+0042369 


-0.97 


G03 


17890048 


2007-0+29 


1.94 


Disk 


106 


102 


18292927+0018000 


-1.24 


G03 


17889792 


2007-04^30 


1.15 


Disk 


107 




18293056+0033377 


-0.60 


G03 


17889536 


2007-04^30 


1.46 


Disk 


108 


105 


18293254-0013233 


-2.58 


G03 


17884672 


2007-04-28 


1.36 


BG star 


109 


106 


18293300+0040087 


-1.35 


G03 


17888512 


2007-05-05 


2.16 


Disk 


110 


107 


18293319+0012122 


-1.81 


G03 


17883648 


2007-04-24 


0.34 


BG star 


111 




18293337+0050136 


-0.26 


G03 


17888768 


2007-05-05 


3.80 


Disk 


112 


109 


18293381+0053118 


-2.06 


G03 


17886720 


2007-0+24 


0.29 


BG star 


113 


110 


18293561+0035038 


-1.53 


G03 


17883904 


2007-04-29 


5.76 


Cold disk 


114 


111 


18293619+0042167 


-0.92 


G03 


17883904 


2007-04-29 


7.52 


Cold disk 


115 


114 


18293672+0047579 


-1.00 


G03 


17886976 


2007-04-29 


3.19 


Disk 


116 


116 


18293882+0044380 


-1.23 


G03 


17890048 


2007-04-29 


1.55 


Disk 


117 


118 


18294020+0015131 


0.87 


G03 


17883392 


2007-04-29 


2.72 


Disk 


118 




18294067-0007033 


-2.53 


G03 


17884672 


2007-04-28 


0.38 


BG star 


119 


119 


18294121+0049020 


-1.39 


G03 


17889280 


2007-04^24 


1.55 


Disk 


120 


123 


18294168+0044270 


-1.42 


G03 


17884416 


2007-04-30 


3.13 


Disk 


121 




18294301-0016083 


-2.18 


G03 


17890304 


2006-10-25 


0.33 


BG star 


122 


126 


18294410+0033561 


-1.68 


G03 


17886464 


2007-04-29 


10.92 


Cold disk 


123 


129 


18294503+0035266 


-1.26 


G03 


17886976 


2007-04-29 


3.71 


Disk 


124 


132 


18294725+0039556 


-1.32 


G03 


17887232 


2007-04-27 


2.22 


Disk 


125 


133 


18294726+0032230 


-1.31 


G03 


17886976 


2007-04^29 


1.50 


Disk 


126 


140 


18294962+0050528 


-2.00 


G03 


17890816 


2007-04^24 


0.35 


BG star 


127 


143 


18295001+0051015 


-1.68 


G03 


17889280 


2007-04-24 


1.23 


Disk 


128 




18295012+0027229 


0.14 


G03 


17888768 


2007-05-05 


2.60 


Galaxy 


129 


144 


18295016+0056081 


-0.76 


G03 


17886208 


2007-09-30 


1.86 


Disk 


130 


145 


18295041+0043437 


-1.32 


G03 


17886720 


2007-04-24 


2.14 


Disk 


131 


148 


18295130+0027479 


-2.42 


G03 


17889024 


2007-04-25 


2.09 


Disk 


132 


149 


18295206+0036436 


-0.06 


G03 


17886464 


2007-04^29 


4.04 


Emb 


133 




18295240+0035527 


-0.33 


G03 


17884416 


2007-04^30 


3.07 


Emb 


134 


153 


18295244+0031496 


-0.59 


G03 


17889792 


2007-04-30 


1.48 


Disk 


135 


154 


18295252+0036116 


0.78 


G03 


17883392 


2007-04-29 


24.39 


Emb 


136 


156 


18295304+0040105 


-1.16 


G03 


17889536 


2007-04-30 


1.50 


Disk 


137 


157 


18295305+0036065 


0.07 


G03 


17882624 


2006-10-23 


10.57 


Disk 


138 


161 


18295322+0033129 


-1.99 


G03 


17886720 


2007-04-24 


2.35 


BG star 


139 


165 


18295422+0045076 


-2.35 


G03 


17889024 


2007-04^25 


1.03 


Disk 
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140 


166 


18295435- 


h0036015 


-0.19 


G03 


17887232 


2007-04-27 


9.65 


Emb 


142 


172 


18295592- 


h0040150 


-0.85 


G03 


17888000 


2007-04-30 


1.21 


Disk 


143 


173 


18295620+0033391 


-0.54 


G03 


17889024 


2007-04-25 


1.21 


Disk 


144 


177 


18295701- 


(-0033001 


-1.15 


G03 


17882880 


2007-04-28 


0.96 


Disk 


145 


178 


18295714- 


h0033185 


-0.82 


G03 


17883648 


2007-04-24 


1.26 


Disk 


146 


182 


18295772- 


h0114057 


-0.44 


c2d 


9828352 


2004-09-02 


1.41 


Disk 


147 


189 


18295872- 


h0036205 


-1.11 


G03 


17888000 


2007-04-30 


0.72 


Disk 


148 


206 


18300178- 


h0032162 


-1.25 


G03 


17886976 


2007-04-29 


1.90 


Disk 


149 


210 


18300350- 


h0023450 


-1.22 


G03 


17888000 


2007-04-30 


1.81 


Disk 



''As in lOIiveira et al.1 l|2009h 
^From lHarvev et al ] l|2007bl ') 

"^^T^ilXr tietween the 2MASS K (2/im) and the MIPSl (24Atm) bands llHarvev et al.ll2007d) 

'*'G03' means the object was observed as part of the Spitzer Space Telescope's cycle 3 program #30223, PI: K. 
Pontoppidan; 'c2d' means the object was observed as part of the c2d 2nd Look program, PI: N. Evans 

Tlux ratio between 30 and 13 fim 



-41 - 



Table 2. Characteristics of YSOs in Serpens 



# 


Emb 


PAH** 


S9.8 


S11.3 


peak 


S19.0 


S23.8 


peak 


Cluster 


Spectral Type 


I 








2.13 


9 


2 00 


1.52 


X . ijyj 




K2 








1 4Q 


J- .0(J 


1 50 


1 33 


1.23 


1.26 




MO 








9 1 9 


1 7*^ 

A., i 


9 1 




1 .00 




±J 


XvO 


7 






-L .00 


1 9^^ 


J. .00 


J. .Ou 


1 '^4 


± .ou 


T) 


ivrn 


o 






1 1 ^ 


1 43 


1.14 


1.21 


1 45 


82 


Y) 




q 






1.64 


1.48 


1 68 

J- ."JO 


2.24 


1 Q7 


2.27 


D 




1 n 






1.25 


1.29 


1.28 


1 47 

J- .4: ( 


1.29 


1 51 

X .(J X 


D 
















1 33 

J- .00 


1 87 


1.29 


D 




14 


- 


- 


1.45 


1.38 


1.45 


1.65 


1.43 


1.53 




M2 


15 


— 


— 


2.67 


2.19 


2.68 


2.59 


2.07 


2.46 


D 
















n qn 


85 


1 1 3 

X . XO 






91 








1 1 7 


1 1 




1 "^7 


1 '^4 

X .O'-t 






94 












-L .U;? 


1 9Q 


X .ou 






9B 


Y 

1 




















97 


Y 




















9R 


Y 




















9Q 








1 i^Q 


J. .Oi7 


-L .(JU 




1 7*^ 
X. 1 




MP 

iVXZi 


ou 








1 '^S 

± . 00 




1 99 
X . zz 


X . xo 


X .00 




Ml 

iVXX 


O -L 






J. . 1 u 


1 79 


1 7fi 

J. . 1 u 


QQ 


Q9 


QQ 
u. y y 






'^9 














X . XtJ 


n 79 






ou 

'^^^ 






1 47 

J- .'± ( 


1 35 

J- .0(J 


1 48 


1.41 


1 31 

X .tJX 


1.41 




K5 


oo 


Y 




















40 






1 f)6 


1.20 


1 OQ 

J- .\Jt7 


1.29 


1.26 


1.25 


E 


M7 


41 






i. (0 


i. ( ( 


1 OQ 


O.40 


/lie 
4.18 


0. /4 






42 


















E 




43 






1.32 


1.40 


1.35 


1.18 


1.18 


1.28 




M0.5 


44 


Y 
















C 


M6 


45 


Y 
















C 




47 


Y 
















E 


M5 


48 






1.22 


1.21 


1.26 


1.12 


1.05 


1.18 


E 


M5.5 


49 


Y 
















C 




50 


Y 
















C 




51 


Y 




















52 




Y 














E 


G3 


53 






1.19 


1.32 


1.25 


1.26 


1.04 


1.18 




M2.5 


54 






1.13 


1.26 


1.21 


1.19 


1.47 


1.07 


C 
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II 

# 


rLmb 


TZ> A xja 




o 

Oil. 3 


qlO/^m 
peak 


O19.0 


'523.8 


q20)Lim 

peak 


duster 


bpectral iype 


55 






1.88 


1.68 


1.89 


1.53 


1.18 


1.35 


E 


K2 


56 






1.56 


1.49 


1.55 


1.32 


1.23 


1.41 






57 






1.39 


1.47 


1.43 


2.01 


1.81 


1.90 


C 




58 






1.09 


1.15 


1.16 


1.13 


1.20 


1.33 


C 


K7 


59 


















c 




60 






1.00 


1.08 


1.01 


1.27 


1.25 


1.26 


c 


MO. 5 


61 






1.24 


1.25 


1.25 


1.28 


1.32 


1.26 


c 


MO 


62 












1.21 


1.31 


1.26 


c 


KO 


63 


Y 
















c 




64 












1.41 


1.27 


1.36 


c 




65 






1.00 


1.10 


1.00 


1.21 


1.12 


1.06 


c 




66 






1.19 


1.18 


1.20 


1.15 


1.07 


1.11 


E 


K5 


67 


Y 
















c 




69 












2.14 


2.17 


2.27 


B 




70 












1.29 


1.36 


1.31 


E 


A3 


71 






1.50 


1.49 


1.48 


2.48 


2.86 


2.51 


E 


M3 


72 


Y 
















c 




73 


Y 
















c 




74 






0.99 


1.24 


1.03 


1.44 


1.33 


1.55 


B 




75 






1.37 


1.37 


1.38 


1.64 


1.35 


1.50 






76 






1.58 


1.51 


1.61 


1.83 


1.78 


1.75 




Ml 


79 


Y 
















c 




80 






1.27 


1.39 


1.29 


0.96 


0.92 


0.98 


B 




81 






1.59 


1.52 


1.66 


1.50 


1.28 


1.42 




M5 








J. .o± 




1 97 


1 14 




1 1 ^ 


E 

Mi 


Mfl 


83 


















E 




86 






1.80 


1.61 


1.83 


1.28 


0.98 


1.23 




M5.5 


87 


















B 


M4 


88 






2.58 


2.75 


2.53 


1.70 


1.41 


1.47 




M0.5 


89 






1.10 


1.21 


1.08 


1.62 


1.68 


1.71 


B 


K5 


90 






1.19 


1.19 


1.21 


1.54 


1.62 


1.67 


B 




91 


Y 
















E 


K7 


92 






1.56 


1.60 


1.63 


1.92 


1.66 


1.97 


E 


MO 


94 






1.14 


0.85 


0.99 


0.46 


0.53 


0.54 






96 






1.35 


1.23 


1.34 


1.53 


1.48 


1.41 


E 


Ml 


97 




Y 
















M2 


98 




Y 
















A3 



-43 - 



Table 2 — Continued 



II 

# 


Jiiinb 


TD A Tja 

rKn 




011.3 


qlOjLim 

peak 


f5l9.0 


c 

023.8 


peak 


Oluster 


opectral iype 


100 






1.39 


1.47 


1.39 


1.90 


1.56 


1.89 






101 






1.60 


1.55 


1.63 


1.95 


1.65 


1.92 


B 




102 


Y 
















B 




103 






1.27 


1.53 


1.31 


1.30 


1.28 


1.16 






104 






1.50 


1.67 


1.53 


1.96 


1.70 


1.90 






105 






0.87 


0.87 


0.80 


3.84 


3.15 


3.30 






106 






1.19 


1.28 


1.19 


1.22 


1.42 


1.61 




M3 


107 






1.07 


1.04 


1.06 


1.24 


1.16 


1.37 






109 












1.58 


1.31 


1.51 


B 




111 






1.24 


1.32 


1.29 


0.61 


1.11 


1.23 






113 






1.10 


1.25 


1.12 


1.76 


1.82 


1.91 




K7 


114 






1.14 


1.21 


1.16 


1.90 


1.99 


2.05 


B 


F9 


115 






1.43 


1.36 


1.45 


2.22 


2.06 


2.11 




MO. 5 


116 












1.56 


1.51 


1.49 


B 




117 






1.37 


1.38 


1.40 


1.87 


1.68 


1.81 




K2 


119 






1.20 


1.23 


1.25 


1.22 


0.82 


1.05 


B 


K7 


120 




Y 








1.70 


1.52 


1.64 


B 


A2 


122 






1.22 


1.17 


1.21 


0.90 


1.08 


1.09 


B 


MO 


123 






1.87 


1.70 


1.90 


1.97 


1.63 


1.95 


B 


MO 


124 






1.09 


1.15 


1.13 


1.06 


1.15 


1.08 


B 


MO 


125 






1.57 


1.50 


1.61 


2.06 


1.86 


2.02 


B 


MO 


127 












1.85 


1.71 


2.00 


B 


M2 


129 






1.19 


1.16 


1.20 


1.11 


0.95 


1.15 






130 












1.39 


1.24 


1.23 




K6 






Y 


















132 


Y 
















B 




133 


Y 
















B 




134 






0.76 


0.88 


0.73 


1.25 


1.33 


0.89 


B 




135 


Y 
















B 




136 






1.08 


0.98 


1.07 


1.30 


1.38 


1.54 


B 




137 






1.27 


1.24 


1.28 


1.99 


2.06 


2.01 


B 




139 












1.88 


1.50 


1.76 




A4 


140 


Y 
















B 




142 






1.50 


1.39 


1.49 


1.26 


1.13 


1.16 


B 


M4 


143 












1.42 


1.26 


1.47 


B 




144 






0.98 


1.11 


0.99 


1.30 


1.23 


1.32 


B 




145 












1.27 


1.26 


1.32 


B 


G2.5 
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Table 2 — Continued 



# Emb 


PAH^ S9.8 


S11.3 


qlO/xm 
peak 


S19.0 


S23.8 


q20/im 
peak 


Cluster 


Spectral Type 


146 


1.23 


1.39 


1.25 


1.28 


1.15 


1.22 


A 


M4 


147 


1.41 


1.51 


1.46 


1.40 


1.23 


1.42 


B 




148 


1.81 


1.67 


1.83 


1.76 


1.51 


1.79 


B 


K7 


149 


1.90 


1.65 


1.89 


1.47 


1.31 


1.43 




MO 



'^All disk sources with PAH in emission show the features at 6.2, 7.7, 8.6, 11.2 and 12.8 fim 



Table 3. Number of Objects in Tight Groups & Associations 





Embedded 


Disks 


Total 


Cluster A 





1 


1^ 


Cluster B 


5 


26 


31 


Group C 


9 


10 


19 


Group D 





8 


8 


Group E 


3 


12 


15 



^The Serpens Core, or Cluster A, is not 
complete in these observations. 



Table 4: Statistics of the YSO populations in Serpens and Taurus 





Embedded 




Disks 








PAH 


both 10 and 20 //m only 20 //m 


featureless 


Serpens 


18.3% 


2.6% 


60.0% 13.9% 


4.3% 


Taurus 


20.1% 


2.9% 


57.6% 3.6% 


18.7% 
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A. Background Sources 

Here the background contamination discussed in §3.11 is presented in 3 categories. 



A.l. Background stars 

The IRS spectra of bright background stars, when seen through a molecular cloud, show 
the silicate feature at 10 /im in absorp tion and not in e missi on as for disk sources (see § l3.3.3p . 
on top of a falling infrared spectrum. lOliveira et al.l (120091 ) identified candidate background 
stars based on optical spectra for 78 objects in this sample. After determining spectral 
types, extinctions and luminosities, those objects were placed in a Hert zsprung-Russell dia- 



gram (HRD) an d compared with 2 sets of isochrones and mass tracks (IBaraffe et al.lll998 



Siess et al.ll2000l ). 20 out of those 7 8 objects proved to be much too luminous to be at the 
distance of Serpens {d = 259 ± 13pc, IStraizys et al.lll996l ). Those objects were then classified 
as background sources, typically asymptotic giant branch ( AGS) stars w i th du sty shells. 
Interestingly, 18 of the 20 objects classified as background in lOliveira et al.l (120091 ) have IRS 
spectra characteristic of background sources; the other two objects, #81 and 86, show sili- 
cates in emission (just like those seen in circumstellar disks) and were then re-classified as 
YSOs. 

Besides these 18 objects, another 8 objects in this sample with falling mid-IR spectra 
show the 10 p^m silicate feature in absorption (objects #19, 25, 37, 78, 95, 126 were not 
observed by lOliveira et al.ll2009l ). leading to their classification as background stars. Alter- 
natively, those 8 objects could be highly extincted YSOs, with sihcate features in absorption 
arising from such extinction. However, using the dense cloud dust properties, as in ^ [37375| 
an extinction of 10 mag would produce a rg.y of only 0.14 (see also lChiar et al.ll2007l ). Since 
these objects do have silicate features in absorption, it is understood that they have high 
extinction (considerably higher than what is found for the cloud at their location). Thus, it 
is unlikely that they are actually highly extincted cloud members. Our working assumption 
is that these are background objects. These 26 spectra are shown in Figure | 



A. 2. Background galaxies 

Five of the sources with clear PAH features show wavelength shifts in their line positions 
compared with galactic sources. These features thus allow us to identify them as background 
galaxies in this sample. PAHs are common features in most nearby galaxies with ongoing or 
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5 10 15 30 35 30 3& 5 10 15 30 S5 30 3& 



Fig. 16. — IRS spectra of the background stars toward Serpens with the characteristic 10 /zm 
sihcate absorption feature and faUing SED. In grey are the original spectra, while in black a 
binned version is overplotted. The object # is indicated at the top right of each panel. 



recent star-f ormation f|Smith et al.l 120071 ). and some contamination by background galaxies 
is expected ( Papovich et al.l 2004 ). 



By identifying PAH features from their shape and pattern (jTielensll2008l and references 
therein), the observed wavelength (Aobs) can be related to the rest-frame wavelength (Arest) 
of each feature through a redshift {z = — 1.0). Between 3 (#12) and 5 (#11 and 68) 
PAH features are identified in each spectrum, and the redshift of each line is calculated. For 



a given object, the z values derived from the various hues agree within 2% and the mean of 
these values is taken as the redshift for that source. These 5 redshifted galaxies can be seen, 
in rest-frame wavelength, in Figure [T71 




Fig. 17. — IRS spectra of the background galaxies toward Serpens, shown in rest-frame 
wavelength. The dashed lines indicate the rest-frame PAH features at 6.2, 7.7, 8.6, 11.2, 
12.8 and 16.4 \xm. The object # and z are indicated in each panel. 

Star-forming galaxies show strong correlations between the extinction-corrected Paschen 
a luminosity coming from HII re gions and the luni inosity at 24 /im due to the emission of 



dust heated by hot, young stars. iRieke et al.l (120091 ) obtained power-law fits to the relation 



between the star formation rate (SFR) and the observed 24 /zm flux densities (/24,06s), pa- 
rameterized with respect to z to take into account the fact that the rest-frame wavelengths 
sampled change with redshift (Eq. 14 in their paper). In order to estimate the infrared lumi- 
nosity (L(TIR)), we use their best-fit relations between the (rest-frame) 24 /xm luminosity, 
L(TIR) and SFR (Eqs. 11 and 25). 
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We have calculated SFRqj and L(TIR) for the 5 background galaxies contaminating this 
sample. These results can be found in Table [51 L(TIR) ranges from a few times 10^^ to a few 
times 10^^ Lq, indicating that these galaxies are typical so-called luminous infrared galaxies 
(LIRGs, L(TIR) > 10" Lq), while two are ultra-luminous (ULIRGs, L(TIR) > lO^^ Lq). 
The latter are a class of heavily obscured g alaxies having enormous instantaneous SFRs that 
are often associated with merging galaxies (jSanders fc Mirabellll996l ). These galaxies are very 
rare locally and become more common at high redshift. Th eir space density rises by about 
a factor of 10 from 2; = to ^ ~ 0.5 (IMagnelli et al.l 120091 ) . consistent with the relatively 
high redshifts found {z ~ 0.3 — 0.5). Furthermore, due to its high SFR an d the erapirica l 
correlation between the infrared and radio emission from starburst galaxies (ICondonlll992l ). 
the ULIRG #12 at z = 0.46 can be identified with a p reviously unidentifie d radio source in 
the NRAO VIA Sky Survey (NVSS J182817+001604. ICondon et al.lll998h . 



The inferred SF Rs range between 11 and 731 M© yr ^, compared to a Milky Way SFR of 

, and referenc es therein). For further reference, 



about 3-5 Mq yr ^ ( jPrantzos fc Aubert 



995 



the SFR of Serpens is 5.7 x 10 ^ Mq yr ^ (lEvans et al.ll2009l ). Even though the star formation 



observed in those galaxies is due to high mass star-forming regions such as Orion, the high 
SFRs in the LIRGs and ULIRGs found are equivalent to 10^ — 10^ molecular clouds such as 
Serpens! 



A. 3. High Ionization Object 

One of the objects, #17 (SSTc2dJ18282720+0044450), in this sample presents very 
strong emission lines that do not seem compatible with YSOs. Figure [TK] shows a blow-up 
of the spectrum, identifying the high ionization lines at 2; = 0. The source itself is highly 
reddened, being barely detected in IRACl (3.5 fim) but very prominent in MIPSl (24 /^m), 
with a 3-24 /xm slope of 1.82. Two types of galactic objects show such high excitation 
lines: dusty planetary nebulae and supernova remnants. The IRS spectrum of this object 



#1 7) resembles those o 



and Ghavamian et a. 



Bernard-Salas et all toO^ . 



shocked ejecta of supernova remnants, as in lSandstrom et al.l (|2009[ ) 



(120091) and that of planetary nebulae such as LMC 78, shown in 



The spectrum also shows PAH emission at 6.2, 7.7, 11.2, and 12.8 /xm. Optical and 
near-IR spectra are necessary to confirm the nature of this object. 



^The SFRs correspond to an IMF according to IChabrieij (|2003h . 
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Fig. 18. — IRS spectra of the planetary nebula candidate(#17) found among this sample, 
presenting high ionization lines. 
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Table 5. Characteristics of the Background Galaxies 



# 


z 


24/im (Jy) 


SFR (Mo yr-i) 


L(TIR) (Lq) 


Type 


11 


0.1051 


0.01800 


16 


1.5x10" 


LIRG 


12^ 


0.4637 


0.00893 


731 


4.8x10^2 


ULIRG 


33 


0.2742 


0.01490 


217 


1.6x1012 


ULIRG 


68 


0.2032 


0.00909 


47 


4.0x10" 


LIRG 


128 


0.1164 


0.01010 


11 


1.1x10" 


LIRG 



'^NVSS J182817+001604 



